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ABSTRACT 

The Kelvin-Helmholtz Instability (KHI) has been observed in the solar atmosphere. Ion-neutral collisions 
may play a relevant role for the growth rate and evolution of the KHI in solar partially ionized plasmas as 
in, e.g., solar prominences. Here, we investigate the linear phase of the KHI at an interface between two 
partially ionized magnetized plasmas in the presence of a shear flow. The effects of ion-neutral collisions and 
compressibility are included in the analysis. We obtain the dispersion relation of the linear modes and perform 
parametric studies of the unstable solutions. We find that in the incompressible case the KHI is present for any 
velocity shear regardless the value of the collision frequency. In the compressible case, the domain of instability 
depends strongly on the plasma parameters, specially the collision frequency and the density contrast. For high 
collision frequencies and low density contrasts the KHI is present for super- Alfvenic velocity shear only. For 
high density contrasts the threshold velocity shear can be reduced to sub-Alfvenic values. For the particular 
case of turbulent plumes in prominences, we conclude that sub-Alfvenic flow velocities can trigger the KHI 
thanks to the ion-neutral coupling. 

Subject headings: Instabilities — Sun: filaments, prominences — Sun: corona — Sun: atmosphere — Magne- 
tohydrodynamics (MHD) 



1. INTRODUCTION 

The Kelvin-Helmholtz Instability (KHI) is a well-known 
magnetohydrodynamic instability that arises at the inter- 
face between two fl uids in relative motio n (see the clas- 
sical textbooks by IChan drasekhad Il961t iDrazin & ReidI 
119811) . The KHI is believed to operate on many astro- 
physi cal plasmas as, e.g., the magneto pause (e. g. Hasegaw3 
1975b , planetary magnetospheres (e.g., lNaganolQ979; Miural 
1984), Eart h's aurora (e.g., iF arrugia et al. 1994), cometary 
tails (e.g., Ershkovich et al.' n986l) , protoplanetary disks 
(e.g., ^G omez & Ostriker 20()5|), jets and outflows (e.g., 
IKeppens et al.1 11999t IBatv & KeppensI |2006|) , among many 
other situations. 

The KHI is also at work in the so lar at mosphere. Recent ob- 
servations by iFouUon et alj (1201 ll) and lOfman & Thomps"onl 
(|201 1 ) confirmed the presence of the KHI in the solar corona. 
Although the KHI has been observed for the first time in 
the corona recently, there is a large number of theoretical 
papers on the KHI in solar plasmas in the literature (e.g. . 



two incompressible plasmas is always unstable in the absence 
of a magnetic field. The presence of a magnetic field com- 
ponent along the flow direction suppresses the KHI for sub- 
Alfvenic velocity sh ear in fully ionized plasmas (see, e.g., 
IChandrasekhaij[l96ll) . However, the situation is more com- 
plicated if the plasma is partially ionized as happens, e.g., in 
prominences. In a partially ionized plasma composed of ions 
and neutrals, the two species behave very differently when a 
magnetic field is applied. Neutrals are insensitive to the mag- 
netic field so that their natural tendency is to be unstable for 
any velocity shear, while ions feel the stabilizing presence of 
the magnetic field. The coupling between neutrals and ions 
through collisions is therefore crucial for the growth rate of 
the KHI and its evolution (see a discussion on the relevance 
of ion -neutral interaction for plasma instabilities in Lehner^ 
[19781) . 

The KHI in partially ionized incompressible plas- 
mas has been studied under different contexts (e.g.. 



2001 



2010; 



Rp 119831; lAndries et al." 120001; 'Andries & Goo ssensI 
Holzwarth et al. 2007; Terradas et al. 2008; Soler et al.l 
Zaqarashvili et al. 2010; Diaz et al. 201 1). More relevant for 
the present investigation are the observations by Berger et al. 
(ir008l 120 ig 12011) and .Ryutova et all (12010) of turbulent 
flows and instabilities in solar prominences, which have been 
interpreted in terms of the Rayleigh-Taylor instability and the 
KHI. The peculiar properties of prominences, i.e., relatively 
cool and dense plasma condensations supported in the corona 
by the magnetic field, make them a very interesting subject 
for the investigation of the KHI. In particular, the fact that the 
prominence plasma is only partially ionized can have impor- 
tant effects on the behavior and evolution of the KHI. 
It is known that a velocity shear at the interface between 

lroberto.soler@wis.kuleuven.be I 



Chhailani & Vvas 1990; Birk & Wiechen 2002; Watson et al 
2004; Shadmehri & Downes 2007, 2008; Kunz 200 ' 
Jones & Downes 201 1). The paper that we take as a referenc e 



for the present investigation is that by Watson et al.l (120041) . 
These authors studied the KHI due to shear flow at the inter- 
face between two incompressible and partially ionized plas- 
mas. They applied their results to weakly ionized outflows 
from massive stars interacting with the ambient interstellar 
material. They concluded that ion-neutral collisions cannot 
suppress the instability of neutrals for sub-Alfvenic velocity 
shear. Thus, a partially ionized incompressible plasma is un- 
st able for any velocity shear. 

IWatson et al.l (12004 ) did not take the effect of compressibil- 
ity into account. In fully ionized plasmas, compressibility can 
have either a stabilizing or destabilizing role depending on 
the density contrast and the orientation of the flow with re- 
spect to the magnetic field direction. For a flow parallel to 
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the magnetic field, compressibility has a stabilizing effect for 
small density contrast and large velocity shear, and a destabi- 
lizing effect for large density contrast and small velocity shear 
(see e.g. iFeier 1964; Sen 1964; Gerwin 1968; Miura 1982; 
lRaelll983l) . In partially ionized plasmas, the role of compress- 
ibility for the KHI was studied by Prialniket al. (1986) and 
lErshkovich et al.l (119861) in the context of comet ionopauses. 
These authors took the flow of neutrals across the interface 
and did not study the case in which the flow of neutrals is 
along the magnetic field direction. To our knowledge, the 
joint effect of compressibility and ion-neutral collisions on the 
KHI has not been studied when the flow of both ions and neu- 
trals is along the magnetic field. 

In this paper we investigate the linear phase of the KHI at 
a sharp interface between two partially ionized plasmas in the 
presence of a shear flow. The flow direction of both ions and 
neutrals is along the magnetic field. Our study focuses on the 
effects of compressibility and ion-neutral collisions. In a par- 
ticular application, we consider typical values of prominence 
plasma parameters. 

This paper is organized as follows. Section |2] contains 
the description of the equilibrium configuration and the ba- 
sic equations. In Section[3]we study the KHI at an interface 
between two plasmas in the coUisionless case, i.e., ion-neutral 
collisions are neglected. Later, we investigate in Section|4]the 
effect of ion-neutral collisions on the growth rates. Section |5] 
contains an application to solar prominence plasmas. Finally, 
our conclusions are given in Section|6] 

2. EQUILIBRIUM AND BASIC EQUATIONS 
2.1. Equilibrium configuration 

Our equilibrium is composed of two partially ionized and 
magnetized homogeneous plasmas separated by a sharp in- 
terface. We consider a hydrogen plasma composed of ions 
(protons), electrons, and neutral atoms. We use Cartesian co- 
ordinates. The x-direction is normal to the interface, so that 
the plane x = coincides with the interface. We denote by the 
indices 1 and 2 the regions corresponding to x < and x > 0, 
respectively. In the two regions, the equilibrium magnetic 
field is straight and along the z-direction, namely B - Be., 
with B constant. We denote by Bi and B2 the magnetic field 
strength in the two regions. The equilibrium plasma param- 
eters are homogeneous and constant in both regions and the 
condition of pressure balance is satisfied at the interface. 

In the equilibrium the plasma species flow along the mag- 
netic field with a constant flow velocity. We consider that all 
the species in the plasma flow with the same velocity. The 
equilibrium flow is Vi = Vie, and V2 = ¥267. We denote by 
Vi and V2 the flow velocities in regions 1 and 2, respectively. 
The flow velocities can be different on both sides of the inter- 
face, i.e., Vi + Vo- Hence, there is a velocity shear that can 
trigger the KHI. 

2.2. Governing equations 

Plasma dynamics are governed by the momentum equa- 
tion, the energy equation, and the mass conservation equa- 
tion of ions, electrons, and neutrals. These equations con- 
tain terms accounting for collisions between the species. In 
the two-fluid trea t ment of partially ionized p lasmas (see, e.g., 
iBraginskiH 119651; IZaqarashvih et al.l 1201 lah . ions and elec- 
trons are considered together as an ion-electron fluid, while 
neutrals form another fluid that interacts with the ion-electron 
fluid by means of collisions. The generalized Ohm's Law 



along with the induction equation governing the evolution of 
the magnetic field are obtained by neglecting the electron in- 
ertia. In the remaining of this paper the indices i, e, and n refer 
to ions, electrons, and neutrals, respectively. 

In this investigation we restrict ourselves to the study of 
linear perturbations superimposed on the equilibrium state. 
Hence, the governing equations are linearized. The set of 
coupled differential equations governing li near perturbations 
from t he equilibrium state are (see, e.g., IZaqarashvili et alj 

lIoTTah 

Pi f $ + V ■ VVi] = - Vpie + - (V X b) X B 



- Q-in (Vi - V„) , (1) 

PnJ^ +V-VVnj = -Vpn-ai„(v„-Vi), (2) 

^- Vx(Vxb) = Vx(viXB), (3) 
of 

% + V ■ V/,ie = -rAeV ■ V„ (4) 

at 

% + V ■ Vp„ = -rP„V ■ v„, (5) 



where V is the equilibrium flow velocity, Vi = (vi.f, Vi,,, Vi,) and 
Vn = (vn.i, Vay, v^z) are the components of the velocity pertur- 
bation of ions and neutrals, respectively, pi^ and are the 
pressure perturbations of the ion-electron and neutral fluids, 
respectively, b = {bx, by, b^) are the components of the mag- 
netic field perturbation, p-, and p„ are the equilibrium densities 
of ions and neutrals, respectively. Pie and P„ are the equilib- 
rium gas pressure of the ion-electron and neutral fluids, re- 
spectively, is the ion-neutral friction coefficient, yU is the 
magnetic permittivity, and y is the adiabatic index. 

Equations ^ and (|2]l are the linearized momentum equa- 
tions of the ion-electron fluid and neutrals, respectively. Equa- 
tion (O is the linearized induction equation, and Equations (|4|i 
and (|5]i are the linearized energy equations of the ion-electron 
fluid and neutrals, respectively. In this equations we have per- 
formed some simplifications. Collisions of electrons with ions 
and neutrals are neglected because of the small momentum of 
electrons. Adiabatic perturbations are assumed and nonadia- 
batic mechanisms are omitted in Equations (|4| and (|5]). In ad- 
dition, magnetic diffusion terms are neglected in Equation (O. 
These simplifications enable us to study the KHI analytically. 
The effect of the omitted mechanisms may be investigated in 
future works by means of numerical methods. 

We express the ion-neutral friction coefficient, ain, as 

ffin = PiPiiTin, (6) 

where is the ion-neutral colUsion rate coefficient per unit 
mass. The friction coefficient vanishes in both the fully ion- 
ized and fully neutral cases. Instead of using in the re- 
maining of this paper we use the collision frequency, which 
has a more practical physical meaning. Thus we define the 
ion-neutral, Vin, and neutral-ion, Vni, collision frequencies as 

Vin = PiTin, Vni = PnTin- (7) 

Since Vin and Vni are related by p„Vi„ - piV^i, we use 
Vin jp all the following expressions for simplicity (see 
IZaqarashvili et al.ll201 Ibi for a discussion on the meaning of 
the different collision frequencies). For our following analy- 
sis it is useful to define the total density as p = pi -i- pn, where 
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we neglect the contribution of electrons. We define the neu- 
tral, ^n, and ion, ^i, fractions as = Pnip and = Pi/p, with 
+ - 1- The parameter is used here to indicate the 
plasma ionization degree. This parameter ranges from = 
for a fully ionized plasma to - 1 for a neutral gas. 

2.3. Normal mode analysis 

We adopt an approach based on normal modes. We put 
the temporal dependence of the perturbations proportional to 
exp (-i(L)t), with oj the frequency. As the equilibrium is ho- 
mogeneous in the y- and z-directions, we can write all per- 
turbations proportional to exp (^ikyy + ik^zj, with ky and L 
the wavenumbers in y- and z-directions, respectively. Due 
to the effect of the flow, the frequency is shifted so that 
Q = w - k.V is the Doppler-shifted frequency. The La- 
gran gian displaceme nts of ions, ^i, and neutrals, are (see, 
e.g.. lGoossens et aLi ri992) 



(8) 



The condition of continuity of and at x = imposes 
that the values of ky and k^ are the same on both sides of the 
interface. 

Since we consider compressible perturbations, we use as 
our variables Ai = V ■ and An = V ■ i.e., the divergence 
of the ion and neutral displacements, respectively. Then, we 
can combine Equations ([T|i-(|5]i and arrive at the two following 
coupled equations. 



-k],M^i'V' 



PiVin 



yPi^ + B^-ln \ 5x2 



■ki^K^iO. 



PnVjn 



W\t\\x - Pn/Pi and k^ k\^^, k^, and 'V^ given by 



"-±11' 

Q2 



Q (Q -1- ixvin) 



i-2 - J.2 , 1.2 _ 



-Ai 



Tie 



(9) 



(10) 



(11) 



(12) 



be related to the wavelength of the plasma perturbations so 
that L. may be estimated from observations. Alternatively, 
if observations of wavelengths are not available, L, may be 
defined using, e.g., the typical scale of the variation of the 
plasma parameters along the z-direction. Thus, the dimen- 
sionless KHI growth rate, yxHi, is computed as 



Tkhi 



(16) 



where is the imaginary part of the frequency of the unstable 
solution. The dimensionless velocity shear at the interface, 
Av, is 

IVi - V2\ 

Av^— -. (17) 

^min 

Unless otherwise stated, we use this nondimensionalization in 
the remaining of this paper 

3. COLLISIONLESS PLASMA 

First, we investigate the KHI when collisions between ions 
and neutrals are neglected. This preliminary study will enable 
us to understand better the results when the coupling between 
ions and neutrals through collisions is included. By setting 
Vin = 0, Equations (|9]) and ( fTOb become. 



5^ An ,2 



5x2 
q2/ 



5x2 



-^l„An=0., 



with ^2 ,^ and k]_^ now given by 



''in — % ^ 



(^24-Q2)(,2,2^-Q2) 

klcl^ - Q2 



(18) 
(19) 

(20) 
(21) 



In the absence of ion-neutral collisions Equation ( fTsT i governs 
the ion-electron fluid and Equation ( fT9] l governs neutrals. We 
study these two equations separately. 



--kl + li 



Q2 



f.2^2 2 



2^2 



kic. 



--X 



Q2 - klc\^ Q -H /Vi, 



~2 



f2 + /Vir 



Ai 



~2 _ : 

Crp- C' 



Q. + iv\. 



Tie 



Q. + i{l +x)vi 



(13) 



(14) 



-(15) 



where c^j = B^/MPi is the square of the Alfven velocity of 
ions, c2^ = yPielpi and cl^ - yPn/pn are the square of the 
sound velocities of the ion-electron and neutral fluids, respec- 
tively, and c^jg = '-Ai'-sie/'-'-Ai ''sie-' squarc of the cusp 
(or tube) velocity of the ion-electron fluid. The coupled Equa- 
tions (|9]l and ( fTOl l are the basic equations of the present inves- 
tigation. In these equations, subscripts 1 or 2 must be used as 
appropriate when the equations are applied in the two regions 
of the equilibrium. 

In the following sections we express the results in dimen- 
sionless units. We express velocities in units of c^i„ - 
min(csni,Csn2) and distances in units of the typical length- 
scale in the z-direction, L,. In particular appHcations may 



3.1. Instability of neutrals 
The solutions to Equation ( fT9] l for perturbations vanishing 



at X ■ 



boo are 



. ^ I Ani exp(k_Lnix) if X < 0, 
" 1 An2 exp(-fcin2-^) if X > 0, 



(22) 



where Ani and An2 are constants. To obtain the dispersion 
relation we impose the continuity of the normal Lagrangian 
displacement and the neutral pressure perturbation at x = 0. 
The dispersion relation is 



PnlQ^ +Pn2n2^, = 0, 



(23) 



with Tn - kj_ni/k^„2 the compressibility factor of neutrals. 

Let us consider first the incompressible limit, i.e., Csni — > oo 
and Csn2 — > oo so that — > 1. The analytical solution of 
Equation ( |23] ) in the incompressible case is 



PnlVi -l-pn2V2, 

a> — K, 



± ik, Wi - V2\ 



VPnlPn2 



(24) 



Pnl + Pn2 Pnl + Pn2 

Note that Equation ( |24] | is independent of ky. The first term 
on the right-hand side of Equation ( l24b is the real part of co. 
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which corresponds to the Doppler shift due to the flow. The 
second term is the imaginary part of u) and informs us about 
the stability of the solution. The solution with the + sign is 
a growing, unstable mode, while the solution with the - sign 
is a damped disturbance. The unstable solution is present for 
any value of the velocity shear This i s the classical incom- 
pressible hydrodynam ic KHI (see, e.g. IChandrasekhadlT96lt 
iDrazin & ReidlfTOSTl) . 

Now we study the compressible case. We solve Equa- 
tion ( [23T l numerically using dimensionless units. First we con- 
sider the case of low density contrast. We fix p^ilPni - 2 
and compute the growth rate as a function of Av for different 
values of kyjk, (see Figure [Tf a)). We compare these results 
with that in the incompressible limit given by Equation ( |24] | 
(dashed line in Figure [TJ a)). The ratio kyjk, controls the ef- 
fect of compressibility. For kyjk, — > oo both compressible 
and incompressible results agree. This means that compress- 
ibility plays no role when the wave vector in perpendicular 
to the flow direction. For finite kyjk^ compressibility is able 
to suppress the KHI for velocity shear larger than a critical 
value. Compressibility cannot completely suppress the KHI 
for small velocity shear. In the particular case ky - the crit- 
ical velocity shear remains supersonic. 

Next we set k^jk, - 10, so that compressibility has a small 
effect. In Figure [Tib) we display -/khi versus Av for different 
values of the density contrast, PnilPai- The critical velocity 
shear for stabilization gets reduced as the density contrast in- 
creases. We find that for small shear and large density contrast 
the growth rate in the compressible case is larger than the in- 
compressible value. Hence, when the effect of compressibil- 
ity is small and the density contrast is large, compressibility 
destabilizes for small velocity shear and stabilizes for large 
velocity shear These results agree with early investigations 
on the compressible hydrodynamic KHI (see, e.g.. lFejerlll964t 
ISea.l964: .Gerwin 1968; Miura 1982). 

3.2. Instability of the ion- electron fluid 

Now we study the instability of the ion-electron fluid. The 
solutions to Equation (fTST l for perturbations vanishing at x — » 

+00 are 

^ _ I Ail exp(fciiix) if x<Q, 
' ~ \ Ai2 exp(-fcj.i2x) if X > 0, 



(a) 5 





Figure 1. Results for collisionless neutrals. Dimensionless growth rate, 
yKHl: as a function of the dimensionless velocity shear at the interface, Av. 
(a) Results for p^i/pnl = 2 and different values of ky/k-. (b) Results for 
ky/k- = 10 and different values of Pni/Pni- The dashed lines in both panels 
correspond to the incompressible case. In all cases k-L- = n. 
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S k^/k^ = 0.5 \ 







(25) 
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where An and Ai2 are constants. As before, we impose the 
continuity of the normal Lagrangian displacement and the to- 
tal (gas plus magnetic) pressure perturbation at x = 0. The 
dispersion relation is 



(26) 



with Ti - kj_ii/k^i2 the compressibility factor of the ion- 
electron fluid. 

We repeat the same process as for neutrals and study the 
incompressible case first, i.e., Ti ^ I- The solution to Equa- 
tion ( l26l l in the incompressible case is 







P,2/Pn=2 ' 


/ P,2/P„ = '0 ■ 


20 
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O) = ■ 



PnVi +p\2V2, 



Pi I +Pi2 



Bi 



(Pil + Pi2) 



- (Vi - V2Y 



PilPi2 



(Pil +Pi2) 



1/2 



(27) 



Note that Equation dZTl i is independent of ky. In the absence 
of flows. Equation (|27] | reduces to the frequency of the incom- 



Figure 2. Results for the collisionless ion-electron fluid. Dimensionless 
growth rate, 7khI: as a function of the dimensionless velocity shear' at the 
interface, Av. (a) Results for palpa = 2 and different values of kyjk-. (b) 
Results for kyjk^ = 5 and different values of pi2/pii. The dashed lines in 
both panels correspond to the incompressible case. An all computations, 
CAi/Csie = 5 and k-L = n. 
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pressible surface mode driven by magnetic tension, namely 



tij — +k- 



(28) 



In Equation dZTl l the second term on the right-hand side in- 
forms us about the stability of the solution. If the terms within 
the brackets become negative, the solution with the + sign is a 
unstable mode, while the solution with the - sign is a damped 
mode. This is the cla ssical incompressible magnetohydro- 
dynamic KHI (see, e.g. ' Chandrasekharl[T96Tt iDrazin & ReidI 
[198 L) . The threshold velocity shear that triggers the KHI is 



Pnp\2 



(29) 



In the absence of magnetic field, Bi - B2 - and the so- 
lution is unstable for any velocity shear Thus, we recover 
the well-known result that the longitudinal component of the 
magnetic field has a stabilizing effect on the KHI (see, e.g., 
IChandrasekhar 1961). 

Now we investigate the role of compressibility. For sim- 
plicity we have restricted ourselves to the case ca > Cj. First 
we fix the density contrast to pulpn - 2 and explore the so- 
lutions of Equation (|26] | for different values of the ratio kyjk- 
(see Figure |2a)). Incompressible and compressible results 
agree when the ratio k^ jk- is large. As the ratio k^ jk^ de- 
creases, perturbations are more affected by the stabilizing ef- 
fect of compressibility. Compressibility can suppress the KHI 
for velocity shear larger than a critical value. Thus, the solu- 
tion is unstable in a range of velocity shears in between the 
threshold value due to the magnetic field (Equation ( |29] l) and 
the critical value due to compressibility. For ky - the com- 
bined effect of compressibility and the magnetic field com- 
pletely suppresses the KHI. 

Again, we study the case in which compressibility has a 
small effect and the density contrast is large. To do so, we set 
ky/k, = 5 and compute the growth rate versus Av for different 
values of pii/pii (see Figure|2b)). When the density contrast 
is large compressibility destabilizes for small velocity shear 
and stabilizes for large velocity shear in comparison with the 
incompressible limit. The behavior of the ion-electron fluid is 
therefore the same as that of neutrals regarding compressibil- 
ity. 

4. COLLISIONAL PLASMA 

Here we consider the coupling between ions and neutrals 
through collisions. Our purpose is to assess the effect of 
ion-neutral collisions on the results of the collisionless case 
(Section [3]l- We express the collision frequency, Vin, in units 
of the surface wave frequency, defined in Equation (l28T l. 
Thus, the ratio Vin/wk controls the relevance of ion-neutral 
collisions. The larger vin/wk, the more important the effect of 
collisions. 

To find the solutions of the coupled Equations (|9]l and (fTOl l, 
we use Equation ( fTOl i to express Ai in terms of A,, and its 
derivatives as 

Next we use this last expression in Equation (|9]l to obtain an 



equation for An only, namely 



_ (k^ + j-2 - feM + (k^ k^ - k^k^ ]a -0 



with k^ defined as 



r/^„(rAe + BV/^)■ 



(3l) 



(32) 



Equation dSTT l is a linear differential equation with constant 
coefficients. The solution to Equation ( ISTT l is of the form. 



An - A exp (kx) , 



(33) 



with A an arbitrary constant and k a wavenumber to be deter- 
mined. By substituting this expression in Equation ( |3TI ) we 
obtain an algebraic equation for k^, whose solution is 



7 V 



+ k^ 



1/2 
(34) 



Hence, two independent wavenumbers, namely k\ and k\ are 
possible. The general expression of An must contain the con- 
tributions of both k\ and k^ . 

To understand why there are two different values of k^ it 
is instructive assuming a weak coupling between the species, 
so that the quadratic terms in Vin can be neglected and k% x 
0. Then the two independent wavenumbers in Equation yUl 
simplify to 

k+ ~ ^±ie' ^- ~ kj_n- (35) 

The two values of k^ reduce to the wavenumbers of the ion- 
electron fluid and neutrals in the collisionless case, respec- 
tively (see Equations (fTTl i and (fT2l) ). In the general case, 
k^ and these two wavenumbers are coupled due to ion- 
neutral collisions. 

Taking into account the above expressions and the condi- 
tion that perturbations must vanish at jc ^ +00, we write the 
general form for An on both sides of the interface as 

A _ j A^ exp{k+ix) + A^ exp{k-ix) if x < 0, ,,,, 
^"'{A^ exp{-k^2x) + A- exp(-k^2x) if x > 0, 

where A]", A^, and A^ are constants. We need to impose 
appropriate boundary conditions at jc = to obtain the dis- 
persion relation. The boundary conditions are the continuity 
of the normal Lagrangian displacement of ions and neutrals, 
the ion-electron total pressure perturbation, and the neutral to- 
tal pressure perturbation. After expressing these quantities as 
functions of An and its derivatives and applying the bound- 
ary conditions, we get an algebraic system of four equations 
for the constants A|, A^ A^, and Aj . The nontrivial solution 
of the system provides us with the dispersion relation, i.e., 
the determinant formed by the coefficients of the equations 
set equal to zero. The expression of the dispersion relation is 
given in the AppendixlAl 

4. 1 . Incompressible case 

Before embarking on the investigation of the unstable so- 
lutions in the compressible case, we first discuss the incom- 
pressible limit. This is the case studied by , Watson et afl 
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Figure 3. Results in the incompressible case. Dimensionless growth rate, 
yKHli as a function of the dimensionless velocity shear at the interface, Ai', 
in a coUisional plasma with Vini/w|j = Vi„2/'^k = OA, 1, and 10. In all 
computations pa/Pl = 2, fnl = fn2 = 0.5, and k-L^ = n. 



(120041) . In the incompressible case, the dispersion relation 
is 

{pil [Qi (Qi + !>^lVinl) - fc?C^n] 
+ Pi2 [i22 (^2 + !)r2Vin2) - ^^C^.^]} 
X [palQ.1 (ill + iVinl) + Pn2^2 (^2 + iVm2)] 
+ (Pnl ill Vinl + Pn2il2Vin2)^ = 0. (37) 

Using the same notation as I Watson et alj (120041) . our Equa- 
tion dSTT i reverts to their Equation (26). 

The solutions of Equation ( l37T i are plotted in Figure [3] We 
want to assess how the results in the collisionless case change 
when we increase the collision frequency. For this reason, 
we consider three different values of the ratio Vin/wk. In 
the case Vin/wk ^ 1 (see the solid lines corresponding to 
I'in/wk = 0.1) we essentially recover the results of the col- 
lisionless case. There are two unstable solutions or branches, 
the first one related to neutrals and the second one related 
to the ion-electron fluid. The growth rates of the two solu- 
tions are very similar to those indicated in Section [3] for the 
collisionless case. Ion-neutral collisions become important 
when Vin and Wk are of the same order (see the dotted line 
corresponding to Vin/wk - !)■ Due to collisions, the growth 
rate of the two unstable solutions decrease. For vin/wk » 1 
(see the dashed line corresponding to Vin/wk = 10) there is a 
strong coupling between the two fluids, so that the two unsta- 
ble branches are affected by the behavior of both fluids. The 
growth rate of the first unstable branch is much larger than 
that of the second branch when Vin/wk » 1. 

Importantly, we find that in the incompressible case the first 
unstable branch is unstable for any velocity shear regardless 
the value of Vin/wk- The growth rate for small velocity shear 
continuously decreases when vin/wk increases, but stabiliza- 
tion is never achieved. Figure |4] displays the growth rate of 
the first unstable branch as a function of Vin/i^k for small ve- 
locity shear (Av = 0.2). In Figure |4] the dashed line corre- 
sponds to the incompressible case. As for small shear the 
solution is governed by neutrals, this means that ion-neutral 
collisions are not able to stabilize neutrals even in the limit 
Vm/oJk » 1- Here our results are consistent with those of 
IWatson et al.l (120041) . who found that the growth rate of the 
solution related to neutrals is app roximately propo rtional to 
vr' (see details in the Appendix of I Watson et al.ll2004l) . 

4.2. Compressible case 



Figure 4. 7khi as a function of Vin/wk, in a coUisional plasma with Av = 0.2 

for Vini = Vi„2 = Vin, Pl/Pl = 2, CAi/Cjie = 5, Cjie = Csn, fnl = fn2 = 0.5, and 

k^L = 7T. The meaning of the different lines is indicated within the Figure. 
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Figure 5. Dimensionless growth rate, yKHl. as a function of the dimen- 
sionless velocity shear at the inteii'ace, Ai', in a coUisional plasma with (a) 
nnl/'Jk = Vml/i^k = 0.1, (b) v,„i/wk = Vml/t^k = 1, and (c) V,„i/&(k = 
yml/t^k = 10. The meaning of the different lines is indicated within the pan- 
els. In all computations pa/Pl = 2, CAi/Csie = 5, c^ie = Csn, fni = fn2 = 0.5, 
and k-L. = n. 



Here we discuss how the results in the incompressible limit 
are affected by compressibility. Compressibility was not con- 
sidered by .Watson et al.. (.2004,) . For simplicity, we restrict 
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ourselves to the case ca > Cg. We compare the results when 
we increase Vin/wk. In the panels of Figure |5] we plot the 
results for three different values of Vin / with a density con- 
trast of P2/P1 = 2. Different line styles are used depending 
on the value of ky/k.. For comparison purposes we have also 
plotted using dashed lines the equivalent growth rates of the 
incompressible case. We obtain that for ky/k, s> 1 the growth 
rates of the incompressible limit are recovered. The results 
displayed in Figure|5]show a complicated evolution of the un- 
stable modes when Vin/wk is increased. In the following para- 
graphs we summarize the main results contained in Figure |5] 

When Vin/cok ^ 1 (see Figure [5^ for vin/wk = 0.1) we 
obtain two unstable branches as in the incompressible case. 
Now the two branches are affected by compressibility. Com- 
pressibility causes the solutions to be closer to the stability 
threshold, i.e., their growth rates are smaller than in the in- 
compressible case. However, compressibility does not result 
in the complete stabilization of all solutions. In particular the 
solutions with large ky/k, cannot be completely stabilized by 
compressibility. Only the solutions with small ky/k- can be 
stabilized by compressibility for large enough Av. 

When Vin and Wk are of the same order (see Figure |5j' 
for Vin/<Wk - 1) two unstable branches are present when 
ky/k, » 1 only (see the solid lines in FigurejSj) corresponding 
to ky/k. = 10). However, as ky/k, decreases the second branch 
disappears and only the first branch remains. The behavior of 
the first branch is complex because it has two regions of insta- 
bility. The first one takes place for small velocity shear and is 
not visible at the scale of Figure [Sj) because it has very small 
growth rates. The second region of instability is the solution 
visible in Figure |5]3 and has much larger growth rates. It is 
only present for velocity shear larger than a threshold shear. 
The threshold shear becomes super-Alfvenic when ky/k- de- 
creases. This is an important difference compared to the in- 
compressible case. 

Finally, when vin/wk » 1 (see Figure |5}; for vin/wk = 10) 
only one unstable branch is found. It has two regions of insta- 
bility as before, but again the first one is not visible at the scale 
of the Figure |5};. Eventually, the first unstable region disap- 
pears when Vin/wk is increased to a large enough value. The 
second region of instability remains regardless the value of 
Vin/wk and has a threshold velocity shear. The threshold shear 
is super-Alfvenic and is approximately given by Equation ( |29] l 
if the density of ions is replaced by the sum of the densities of 
ions and neutrals. Thus, in the limit Vin/wk — > 00 we obtain 
the behavior of a plasma where ions, electrons, and neutrals 
behave as a single-fluid. In this limit, the linear growth rate is 
the same as in the fully ionized case but replacing the density 
of ions by the sum of densities of ions and neutrals. 

The results discussed in the previous paragraphs are qual- 
itatively displayed in Figure |6] which corresponds to the in- 
stability diagram Vin/fWk vs. Av when the remaining param- 
eters are kept constant. In the top left part of the Figure 
we can see a small region of instability between two stable 
zones. This is the region of instability with very small growth 
rates that is not visible in Figures |5]3,c. Eventually, for large 
enough Vin/wk this region of instability disappears and the 
two stable zones on both sides merge. For the range of pa- 
rameters considered in Figure |6] the two stable regions merge 
for Vin/wk ~ 70. However we must remark that the critical 
value of Vin/wk depends strongly on the equilibrium parame- 
ters so that it may vary substantially when the parameters are 
changed. Hence we have to be cautious since it is not pos- 
sible to give a general value for the ratio Vin/wk to achieve 




Figure 6. Instability diag ram Vi^/wk vs. Av with P2/P1 — 2, CAi/c^ie — 5, 
Csie = Csni ^ni = fn2 = 0.5, ky/k^ = 2, and k,L = n. Tlie shaded areas 
denote regions of stability, while the white zone denotes instability. The three 
horizontal dashed lines correspond to the values of Vin/aik used in the three 
panels of Figure |5] The vertical dotted Hne corresponds to the value of Ai' 
used in Figure|4] 

stabilzation in all cases. 

In Figure |6] we clearly see the transition between the col- 
lisionless case and the single-fluid case. Let us discuss this 
transition in physical terms. For Vm/co^ <sc 1 there is only 
one stable zone due to the combined effect of the magnetic 
field and compressibility. On the left-hand side of the stable 
zone the instability is governed by neutrals, while on the right- 
hand side the instability is mainly due to the ion-electron fluid. 
Thus, the behavior is the same as in the collisionless case. As 
Vin/wk increases the effect of collisions is, first, to destabilize 
the system so that the stable region is removed. Hence for 
a weak coupling between the species the stabilizing roles of 
compressibility and the magnetic field are overridden. This is 
mainly a consequence of the influence of neutrals, which are 
insentitive to the magnetic field, on ions. As Vi,i/wk keeps in- 
creasing and the coupling between species gets stronger two 
additional stable regions appear at low velocity shear This 
first stable region is due to the fact that ion-neutral collisions 
are able to stabilize neutrals for small velocity shear. The sec- 
ond stable zone appears due to the combined effect of col- 
lisions and the magnetic field on the ion-electron fluid. For 
Vin/wk ^ 1 ion-neutral collisions can eventually suppress the 
unstable zone between the two stable regions so that single 
stable zone is present for low velocity shear as in the single- 
fluid limit. 

Compressibility has a strong impact on the KHI in compari- 
son to the incompressible case, specially when we depart from 
the limit Vm/tOk 1. In general, the effect of compressibilty 
gets stronger as ky/k, decreases. For large ky/k- the effect 
of compressibility is weak and the growth rates tend to those 
in the incompressible case. Compressibility can suppress the 
KHI for small velocity shear so that a threshold velocity shear 
appears. We refer again to Figure H] where this result can be 
seen in more detail. Figure H] displays the growth rate of the 
first unstable branch as a function of Vm/cok for small veloc- 
ity shear (Av = 0.2). The smaller ky/k-, and so the stronger 
the effect of compressibility, the smaller the value of Vin/wk 
needed for stabilization. 

Next we study the effect of increasing the density contrast. 
Pi/ Pi, when the effect of compressibility is small. We set 
ky/k^ - 10 and Vm/cok - 10 so that we focus on the case 
Vin/wk » 1. Figure 13 a) displays yKHi versus Av for differ- 
ent values of pi/p\- At the scale of Figure|7|a), compressible 
and incompressible cases give very similar results. The dif- 
ferences are more important for small velocity shear. This 
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Figure?, (a) 7khi. as a function of Av with Vi„i/fcik = Y[^2l'^k = 10 and 
different values of p2/pi. The remaining parameters are the same as in Fig- 
ure|5] (b) Same as panel (a) but for small Av. The dashed lines in panel (a) 
correspond to the incompressible case. 

can be seen in more detail in Figure |2tb), which focuses on 
small values of Av. The results displayed in Figure |7tb) are 
not visible at the scale of Figure|7ja). For small velocity shear 
the growth rate has a complicated dependence on the density 
contrast. Importantly, we find that when the effect of com- 
pressibility is small and high density contrast the threshold 
shear can be reduced to sub-Alfvenic velocities. 

Although the case ca > Cs is appropriate for solar atmo- 
spheric applications, the case c., > ca also deserves some at- 
tention. We have performed several additional computations 
in the case Cs > ca and have assessed the effect of increasing 
the sound velocity. In summary, we find that when the sound 
velocity gets larger the results tend to those of the incompress- 
ible case. Thus, we need to take small values of the ratio kyjk- 
for compressibility to have an impact on the growth rates. In 
the limit Cs — > oo the incompressible case is recovered and the 
growth rates are independent of kyjk,. 

5. APPLICATIONS TO SOLAR PROMINENCES 

Here we apply our results to the case of flows in solar 
prominences. Prominences are cool and dense plasma con- 
densations supported in the solar corona by magnetic forces. 
Due to the relatively low temperature of the prominence 
material (~ lO** K), the prominence plasma is only par- 
tially ionized. Mass flows prominences have been frequently 
repor t ed (e.g., I Zirker et al.' 1998"; 'Wang* '1999"; 'Kucera et alj 
l2()03l;lLinetalf 2003: Okamoto et al. 2007; Ahn et al. 20lC 
Therefore, prominences are a good subject for an application 
of the theory developed in the present paper 

We perform two different applications. First, we consider 
the case of turbulent upflows in quiescen t pr ominences ob- 
served bv lBerger et alj d2008ll20Toll20nh and llvutova et alj 
( I2OIOI) . Later, we study flows along thin threads of promi- 
nence barbs. Here we use dimensional units. We take and 



L, as the typical length scales for the perturbations en the y 
and z-directions, respectively. We compute ky and k, as 



k - - 



(38) 



The expression for the ion-neu tral collision frequ ency, Vi, 
a hydrogen plasma is (see, e.g.. lSoler et al.ll2009t) 



m 



Pi 
2mp ^ 



\6kYiT 



-O-u 



(39) 



where T is the temperature, is the proton mass, k^ is the 
Boltzmann constant, and crjn is the collisional cross section 
((Tin ~ 5 X 10"''^ m^ for hydrogen). In Equation ( [39l ) a strong 
thermal coupling between the species is assumed so that the 
temperature of both ions and neutrals is the same. Finally, we 
define the dimensional KHI timescale, tkhi, as 



1 



Tkhi 



Tkhi 



(40) 



where now 7khi is the dimensional growth rate. In our anal- 
ysis, we have restricted ourselves to the linear stage of the 
KHI. This stage represents the initial phase in which the un- 
stable mode grows after its excitation. Plasma perturbations 
become important during the subsequent nonlinear evolution 
of the KHI. However, the growth rate of the linear phase pro- 
vides us with a reasonable estimation of the typical timescale 
on which the KHI would be observable. 

5.1. Turbulent plumes 

iBerger et al.1 (l2008l I20T0I I20T1I) and lRvutova eTaH (l20Toh 
reported on the existence of turbulent flows and instabilities 
in quiescent prominences. The observations have been in- 
terpreted in terms of Ray leigh-Taylor and Kelvin-He lmholtz 
instabilities ( Be rger et al.1 120101: iRvutovaetal] 1201 01) . Here 
we are interested in the observations of turbulent flows called 
dark plumes (see a detailed explanation of the phenomenon in 
IBerger et al.ll20T^ 

As reported bv lBerger et al.l (120101) . dark plumes are turbu- 
lent upflows in prominences which usually develop Kelvin- 
Helmholtz vortex rolls. In Ca II H-line, plumes are seen dark 
in contrast to the prominence material, which suggests that 
the plasma in the plumes is hotter and probably less dense 
than the prominence material. The width of the plumes ranges 
between 0.5 Mm to 6 Mm and their maximum heights are be- 
tween 11 Mm and 17 Mm. The mean ascent speed is about 
15 km s"' although velocities up to 30 km s"' are measured. 
The typical plume lifetime is between 400 s and 890 s. 

In the following analysis, medium 1 represents the promi- 
nence plasma and medium 2 the plume. We express all quan- 
tities in dimensional units. We use typical parameters of qui- 
escent prominences: Bi = Z?2 = 5 G, pi = 5 X 10"" kg m"^, 
and Ti = 10"* K. The exact ionization degree of the promi- 
nence plasma is unknown (see recent estiniation s in, e.g., 
iGouttebroze & Labrossd2009HLabrosse et al.l2010 ). We take 
^ni - 0.5, meaning that half of the prominence material is 
neutral. For the plume density, we take a contrast of 10 with 
respect to the prominence plasma, i.e., p2 = pi /lO. The plume 
temperature, T2, is computed from the pressure balance con- 
dition, which gives T2 = 10^' K. At the plume temperature 
we assume the plasma is largely ionized and take = 0.1. 
For simplicity, we assume the prominence plasma is at rest, 
i.e., Vi = 0, and use the plume velocity, V2 = Vpiume, as a 
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Figure 8. (a) Application to turbulent plumes. KHI timescale, tkhi. versus 
plume velocity, Vpiumc- The different lines correspond to different values of 
L- indicated within the Figure, (b) Application to thread in barbs, tkhi versus 
the prominence thread ionization degree, fni, for different values of Vthicad- 
in both cases, the remaining parameters are given in the text. 

parameter. We re late Ly to the plume width. According to the 
observations (see lBerger et al.ll2010l) . the width of the plumes 
is between 0.5 Mm to 6 Mm. For the present application we 
take Ly = 1 Mm. L, is more difficult to estimate from the 
observations, but according to the size of the vortex rolls de- 
veloped in plumes, L~ is probably of the same order or smaller 
than Ly. We use L, as a parameter 

We display in Figure |8la) the dimensional KHI timescale, 
Tkhi, as a function of the plume velocity, ypiume, for different 
values of L-. We obtain that tkhi increases as L- decreases. 
This result is consistent with the idea that the smaller L-, and 
so the larger k,, the stronger the stabilizing effect of the mag- 
netic field. Also, we see that for observed plume velocities, 
i.e., Vpiume ^ 30 km s"', tkhi is longer than the observed 
plume lifetimes. For Vpiume > 40 km s"', approximately, tkhi 
becomes of the same order as the observed plume lifetimes. 
However, flow velocities of 40 km s"' are about two times 
larger than observed plume velocities. 

The fact that we have to take flow velocties about two times 
larger than the observed ones to obtain realistic growth rates 
is quite a reasonable result considering the simplicity of our 
model. We must be aware that the present configuration, i.e., 
a magnetic interface, might miss relevant effects for plume 
dynamics. For example, body modes are not described in our 
model and may be important for the KHI. To take body modes 
into account we should consider slab or cylindrical models. 
This is a task to do in forthcoming works after this first inves- 
tigation. In addition, several effects not considered here may 
reduce the KHI timescale for smaller velocities. For example, 
to take the flow of neutrals not strictly parallel to the magnetic 
field direction may decrease the KHI timescale to values com- 
patible with the observed instability and plume lifetimes (see 
iPrialnik eral...l986;,Ershkovich et al..l986.) . A different mag- 



netic configuration may also influence the KHI growth rates. 
However, the important result we want to stress here is that 
the KHI is present for flow velocities which are of the same 
order as those observed. This is due to the neutral-ion cou- 
pling, since in our configuration a fully ionized plasma would 
be stable for such velocities. 

5.2. Threads in prominence barbs 

High resolution observations show that prominences are 
form by a myriad of thin and long field-aligned ribbons, usu- 
ally called threads. Threads are the basic sub-structures of 
prominences. They outline the magnetic field from the barbs 
to the spines of prominences (Engvold 2008). The observed 
widths of the threads are between 100 km and 600 km, while 
their lengths usually range between 3,000 km and 28,000 km, 
respectively. Mass flows along th reads have been frequently 
reported (see the recent review by lLabrosse e t al. 2010). The 
typical flow velocities are less than 30 km s"' in quiescent 
prominences, although larger values up to 40-50 km s"' have 
been observed in active region prominences. A recent re view 
on the properties of prominence threads can be found in iLinI 

dloTl . 

The lifetime of threads in Wa observations is about 10- 
20 minutes, before th ey eventually disappear (e.g., Lin 20041 
iLineLaLl EOOS", ^2009). The cause of this disappearance is 
unknown, although several explanations invoking instabilities 
have been proposed (e.g., Zaq arashvili et al. 2010; Sole r et alj 
1201 ll) . In this application, medium 1 represents the thread and 
medium 2 the inter-thread plasma. We use B\ - B2 - 5 G, 
pi = 5 X 10"" kg m-^ Ll = 10"* K, p2 = pi/100, and 
^n2 - 0.5. Here we use the prominence ionization degree, 
^ni, as a free parameter The inter-thread temperature, T2, 
is computed from the pressure balance condition. Accord- 
ing to the obs erved widths and lengths of the threads (see, 
e.g.. lLinll20TTh . we take Ly = 200 km and L, = 10,000 km. 
Hence, ky/k^ = 50, meaning that we are in the incompressible 
regime. For simplicity, we assume the inter-thread plasma is 
at rest, i.e., V2 - 0, and use the flow velocity in the thread, 
Vi - Vthread, as a parameter 

The behavior of the KHI timescale, tkhi, with Vthread is sim- 
ilar to that shown in FigureOa) for the case of plumes. How- 
ever, the values of tkhi are significantly larger in the case of 
thread conditions. We display in Figure [8lb) tkhi as a func- 
tion of the thread ionization degree, , for different values 
of Vthread- The ionization degree has little impact on tkhi un- 
less we approach the fully ionized (^ni — > 0) or the neutral 
(^ni — > 1) limits, in which the timescale decreases. The very 
large values of tkhi obtained indicate that, although the KHI 
is present, it can hardly be responsible for the instability of 
threads. This result encourages us to pursue other explana- 
tions in the future. 

6. CONCLUSION 

In this paper we have studied the KHI due to shear flow at 
an interface between two partially ionized plasmas. We have 
focused on the effects of ion-neutral collisions and compress- 
ibility. Our results are summarized as follows. 

In a collisionless plasma, neutrals are unstable for any ve- 
locity shear, while the ion-electron fluid is unstable for super- 
Alfvenic shear only thanks to the stabilizing role of the mag- 
netic field. When ion-neutral collisions are included, ions and 
neutrals get coupled. In the incompressible case, we obtain 
that ion-neutral collisions are not able to stabilize neutrals 
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regardless the value of the collision frequency. The KHl is 
present for any veloc ity shear. Hence, our results confirm the 
previous findings by IWatson et al.l (120041) . In the compress- 
ible case, the domain of instability depends strongly on the 
plasma parameters. In particular, the values of the collision 
frequency and the density contrast are relevant. For high col- 
lision frequencies and low density contrasts neutrals become 
stable and the KHl can be suppressed for sub-Alfvenic veloc- 
ity shear. For high density contrast neutrals are more unstable 
and the KHl can be present for sub-Alfvenic velocity shear. 

We have applied our results to solar prominences. In the 
case of turbulent plumes, we find that the KHl is present 
for flow velocities of the order of those observed. The the- 
oretical KHl timescales are larger than the values consistent 
with the observations, which is probably a consequence of us- 
ing a simplified model. On the contrary, in the case of flow 
along threads we find that the KHl cannot be responsible for 
the thread disappearance, so another mechanism must be in- 
volved. 

Here we have restricted ourselves to the linear phase of 
the KHl. Linear growth rates provide us with characteristic 
timescales for the instability to operate. However, nonlin- 
ear studies are needed to assess the real impact of the insta- 
bility on t he evolution of the plas ma parameters. In a re- 
cent work Jones & DownesI (1201 ll) investigated the nonlin- 
ear evolution of the KHl in a partially ionized plasma in the 
ambipolar-dominated regime. Despite the differences in set- 
up and formalism, it is instructive to compare the results of 
I Jones & Downe s' ("201 1) with our present findings. They used 
a two-dimensional configuration in the xy-plane, with the flow 
and the magnetic field parallel to their y-direction, and the 
interface parallel to the x = plane. No dependence of 
the perturbations in the ir z-direction was considered. There- 
fore, the equilibrium of lJones & Downed (1201 11) is equivalent 
to our model with our k^. - and our L set equal to their 
ky. iJones & DownesI (1201 Ih concluded that the linear growth 



rates in the ambipolar-dominated regime are the same as in 
the ideal case, i.e., in the absence of ambipolar diffusion. In 
our formalism the ambipolar-dominated regime in which ions 
and neutrals are strongly coupled is equivalent to our case 
Vin/wk » 1 . We find that the linear growth rate and the thresh- 
old shear in the limit Vi„/tLtk — » °° (see Section |4~2] | are the 
same as in the fully ionized case but replacing the density of 
ions by the sum of densities of ions and neutrals. Hence, our 
present results and the conclusions of lJones & DownesI (1201 ll) 
seem to be in agreement in this particular limit. Our results 
predict that the growth rate depends on Vin / when we depart 
from the limit Vin/wk — » 

The present investigation may be extended in the future by 
taking into account different flow velocities and directions 
for ions and neutrals, by incorporating the effect of gravity, 
and by including nonideal effects such as nonadiabatic mech- 
anisms or magnetic diffusion. Also, the study of the nonlinear 
evolution of the KHl for arbitrary collision frequency is of 
interest. 
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APPENDIX 



DISPERSION RELATION IN THE COMPRESSIBLE COLLISIONAL CASE 

0, with D (oj) given by the solution of the following 



In the compressible collisional case, the dispersion relation is D (cj) 
determinant. 
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an = C5,i^_i + C3,iC6,i^-i {k\ - ^i„i) , 



Gil 

an 

an 

ai4' 
ail ■■ 
ail' 
a23- 
ai\- 

031 ■■ 



(A2) 
(A3) 
(A4) 

(A5) 
(A6) 
(A7) 
(A8) 
(A9) 

(AlO) 
(All) 
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033 - ^5,2^+2 + C3,2C6,2fc+2 (^+2 ~ ^Inl) ' 
a34 - C$^2k-2 + C3,2C6,2fc-2 {k^2 ~ ^\.nl) ' 

a4i - (Cs.iCvj + C8,i) fe+i + C3,iC5,iC6,ifc+i {k\i - k\^^ , 

a42 = (C5,|C7,l +C8,l)fc-1 + C3,lC5,lC6,|fc-l (fc^l - ^Inl) ' 
^43 - (€5,2^7,2 + C"8,2) k+2 + C3,2C5,2C6,2fc+2 (^+2 ~ ^±02) ' 
fl44 = (C"5,2C7,2 + C8,2) k-2 + C^^2Ci^2C(,2k-2 (k^_2 ~ ^'i.tO) ' 

where parameters C\ - Cg are defined as 

2,2 2 2 

C\ =Pi — 



Q2 + /(l+;^.)v„ 
2^2 ^2 



C2 = 





V 


C3 = 


Qvin 


C4 = 




C5 = 


1 Q 

£12 



' 2 , 2 _ ^.^sie^Ai ^ + iVin \ 
^Ai j^2 Q + j(l+Ar)Vinj' 



2„^2 



J.2 2 ;^2 



C6 = - 
C7 = / 

Cs = - 



1 ^^'(4. + '^' 



Ai) ~ ^' 



2 ~2 
sie Ai 



£1 + j'vii 



£12 



Q + iVin 
^2 



£!(£! + j'Vin) 



(A12) 
(A13) 
(A14) 
(A15) 
(A16) 
(A17) 

(A18) 
(A19) 

(A20) 

(A21) 
(A22) 

(A23) 
(A24) 
(A25) 
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